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spherule bed was originally composed of impact-glass spherules 
ejected from the putative KT impact crater on the Yucatan penin- 
sula. The impact nature of the spherule bed was confirmed by the 
previous identification of shocked quartz and relict impact glass. 
These melt-glass spherules were deposited on a shallow shelf, 
probably located north -northwest of the present sites in the Mendez 
basin [3]. Thick, unstable deposits of these spherules were mobi- 
lized as debris flows at the mouths of canyons in the shelf/slope by 
either tsunamis, earthquakes, or storms (associated with the KT 
impact). 

The debris flows either generated their own turbidity currents 
involving unconsolidated shelf sands, as has been calculated theo- 
retically and observed experimentally [6], or these sands were 
mobilized into turbulent flow simultaneously with the spherules, 
unmixing from the latter by a sieve mechanism operating on rela- 
tive particle sizes and flow densities. An alternative to this single- 
pulse hypothesis is that of two separate episodes, wherein the 
spherule-laden debris flows were separate events that preceded 
the turbidity-flow units slightly in time. Arguing against the lat- 
ter dual-pulse hypothesis is the lack of significant amounts of 
spherule-bearing material incorporated in the base of the sandstone 
unit, which should occur if this unit originated later as a separate 
turbidity flow that scoured the top of the previously deposited debris 
flow. 


deposits. For the reasons stated above, we inteipret the clastic 
sequences at Mimbral and La Lajilla as spherule-bearing, debris-flow 
units overlain by sandstone turbidite units. Furthermore, both of 
these clastic units were mobilized on a distant shelf/slope soon after 
impact by the effects of seismic or oceanic disturbances associated 
with the KT impact event. 
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Either process resulted in debris flows consisting of glass sphe- 
rules that scoured and filled shallow channels in basinal Mendez 
marls, followed closely by sand in turbid suspension that was 
deposited on top of the spherule beds in the neoformed channels. 
This model resembles the modem example of upper Pleistocene 
chaotic silt beds transported by debris flows (analogous to the 
spherule beds) and directly overlain by turbidite sand beds in distal, 
channelized, outer-fan lobes of the Mississippi delta in the Gulf of 
Mexico [7]. As in thi s example, the presence of relatively undeformed 
clasts in the Mexican spherule beds suggests deposition in a non- 
turbulent, debris-flow regime. 

Channel scouring by the debris flow is exemplified by the rip-up 
clasts of underlying Mendez marl contained in this bed. Also, the 
bed is deposited in a typical shallow-channel profile with clearly 
defined edges. Both the spherule unit and the sandstone unit contain 
benthic foraminifers of late Maastrichtian age from a shallow neritic 
environment, whereas benthic taxa are rare in the underlying, deep- 
water, Mendez marl [3,8]. This strongly suggests that these clastic 
units are allochthonous, having been transported from shallow 
water to their present basinal setting. 

The sandstone unit at Mimbral displays sedimentary features 
similar to those that define a turbidite. For comparison, these 
features have been marked as divisions of the Mimbral sandstone 
unit in Fig. la, using the same nomenclature (T^) that subdivides 
a classic turbidite sequence (Fig. lb). 

The final, fine-grained phase of turbidite deposition and re- 
commencement of the hemipelagic rain coincided with the arrival 
of the last portion of the fireball layer at the seafloor, components 
identifying this partial fireball layer are an Ir anomaly and mag- 
nesioferrite spinel crystals [9]. The presence of a capping fireball 
layer limits the time of Formation of this clastic unit to a period of 
up to several months following impact, according to our single-impact 
KT model [9]. 

The almost instantaneous loading of a thick deposit of glass 
spherules (up to 3 mm in diameter) onto a shallow shelf constitutes 
a rather rare event that may, when mobilized, produce unique 


Along with the appreciation of the importance of impact events 
on the Earth’s evolution, there has been increasing speculation that 
energetic collisions have been responsible for processes as varied as 
continental flood basalt eruptions, mantle plumes, continental rift- 
ing, and geomagnetic pole reversals. The link between impacts and 
such geophysical processes was first discussed by Seyfert and Sirkin 

[1] , who suggested that impact- induced mantle plumes could be a 
mechanism for initiating the breakup of plates. Burek and W&nke 

[2] listed correlations between known Cenozoic impacts and geo- 

magnetic field reversals, unconformity ages, shifts in paleotemp- 
eratures, and tectonic episodes. They suggested that major impacts 
could generate shock-induced phase transitions in the upper mantle, 
disrupting a delicately balanced stability down to the core-mantle 
boundary. Rampino and Strothers [3] proposed a quasiperiodic 
correlation between mass extinctions and major continental flood 
basalt volcanism over the last 250 m.y. and attempted to explain it 
in terms of episodic showers of impacting comets. Connections 
between impacts and the internal workings of the Earth are sup- 
ported by correlations of the ages of tektites from strewn fields with 
geomagnetic field reversals [4], and by a reversal associated with 
sediments deposited immediately after the impact that formed the 
Ries Crater [5]. - ~ — - 

A causal link between major impact events and global processes 
would probably require a significant change in the thermal state of , 
* the Earth’s interior, presumably brought about by coupling of im- : 

■ pact energy. One possible mechanism for such energy coupling from 

the surface to the deep interior would be through focusing due to 
axial symmetry. Antipodal focusing of surface and body waves from 
earthquakes is a well-known phenomenon [6], which has previously 
been exploited by seismologists in studies of the Earth’s deep 
; interior [7,8]. Antipodal focusing from impacts on the Moon, Mer- 
| cury, and icy satellites has also been invoked by planetary scientists 
1 to explain unusual surface features opposite some of the large 
: impact structures on these bodies [9,10]. For example, “disrupted 

terrains have been observed antipodal to the Cal oris impact basin on 
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Mercury and Imbrium Basin on the Moon. Very recently there have 
been speculations that antipodal focusing of impact energy withjn 
the mantle may lead to flood basalt and hotspot activity [ 1 1 , 1 2], but 
there has not yet been an attempt at a rigorous model. 

We have proposed a new hypothesis and performed preliminary^ 
proof-of-principle tests for the coupling of energy from major im- 1 
pacts to the mantle by axial focusing of seismic waves. Because of I 
the axial symmetry of the explosive source, the phases and ampli- ; 
tudes are dependent only on ray parameter (or takeoff angle) and are" 
independent of azimuthal angle. For a symmetric and homogeneous: 
Earth, all the seismic energy radiated by the impact at a given takeoff 
angle will be refocused (minus attenuation) on the axis of symmetry,' 
regardless of the number of reflections and refractions it has expe- 
rienced. Mantle material near the axis of symmetry will experience 
more strain cycles with much greater amplitude than elsewhere and^ 
will therefore experience more irreversible heating. The situation is' 
very different than for a giant earthquake, which in addition to; 
having less energy, has an asymmetric focal mechanism and a larger; 
area. It should be noted that our hypothesis is fundamentally differ- _ 
entthan those proposed by many others [e.g., 13-16], which involve ; 
melting and excavation at the impact location. Problems with mod- 
els of this type have been pointed out by Loper [17]. 

We are using two independent proof-of-principle approaches. : 
The first makes use of seismic simulations, which are being per- I 
formed with a realistic Earth model to determine the degree of 
focusing along the axis and to estimate the volume of material, if 
any, that experiences significant irreversible heating. The second 
involves two-dimensional hydrodynamic code simulations to deter- * 
mine the stress history, internal energy, and temperature rise as a 
function of radius along the axis, ----- 

For the preliminary seismic modeling the impact was repre- 
sented as a vertical point force applied at the Earth’s surface as a 
delta function in time. The impactor was assumed to yield an 
impulse of approximately 3 x I0 24 dyne sec, as an estimate of a KT- 
sized impact. Synthetic displacement and strain records were gen- 
erated for such a source by summing the normal modes of the elastic 
Earth model 1066A of Gilbert and Dziewonski [18], We used the 
attenuation profile of the model PREM [19] to determine Q values 
for the 1 066A modes. We assumed a vertically directed point source 
so we included only spheroidal modes in the synthetic seismogram 
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Fig. I. Synthetic strain historiesai and near impact antipode, plotted on same 
scale. 


calculations. Toroidal modes are not excited by a vertical point 
force. All spheroidal modes with periods greater that 45 s were 
summed for the synthetics. 

The synthetic signals yield estimates of the peak strains at any 
location in or on the Earth. Figure 1 shows the difference in strain 
histories on the Earth’s surface at two different angular distances 
from the source. This figure demonstrates the gross effects of antipo- 
dal focusing. Strains at the surface near the antipode (angular 
distance * 180°) are orders of magnitude higher than those over most 
of the rest of the Earth’s surface. Figure 2 plots the peak strains as 
a function of depth beneath the antipode of the impact, from the 
surface to the core-mantle boundary. It can be seen that peak strain 
amplitude varies by 2 orders of magnitude and is largest at the top. 
This would imply that most of the energy is focused at shallow 
depths; however, it may be an artifact of ignoring shorter-period 
modes. This issue will be resolved by adding more modes, and by 
doing body wave calculations. Figure 3 demonstrates that focused 
arrivals have much greater amplitudes than direct arrivals at the 
core-mantle boundary. At the beginning of Fig. 3a, the direct arrival 
is shown for a point directly beneath the impact, followed by a long 
high-amplitude focused trace. For comparison, the strain history at 
the same level on the antipodal axis is shown in Fig. 3b. It can be 
seen to have a larger peak amplitude. 

In addition to seismic modeling, we are in the process of perform- 
ing detailed calculations of a hypothetical impact and subsequent 
wave propagation through the Earth with the CTH [20] strong-shock 
hydrodynamics code. CTH is a one-, two-, and three-dimensional 
multimaterial elastic-plastic Eulerian hydrodynamics code. CTH 
accurately computes the very large material deformations, material 
ejection, and strong-shock generation created by the initial KT 
impact. For these calculations we model the Earth as a central, 
gravitationally stable, stratified body, with a variety of interior 
structure, including inner and outer core and mantle. The materials 
are treated as elastic-plastic. Also, the equations of state associated 



Fig. 2. Peak strain vs. depth along antipodal axis. 
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A = 0°. Depth = 2887 km (above CMB) 
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Fig. 3. Synthetic strain histories just above the core-mantle boundary (CMB) 
on the impact axis(a) beneath the impact and(b) antipodal to the impact. 


with CTH allow thermodynamically consistent phase transitions, so 
that vaporization at the initial impact site can be properly modeled, 
as well as shock-wave modification in the Earth’s interior due to 
solid-solid phase transitions. Alternatively, we have also modeled 
the initial impact site as a general energy source, similar in strategy 
to the seismic model and to the work of Watts et al. [10]. 

CTH can be used to study the long-term evolution of the resulting 
wave structures through the interior of the Earth, with some caveats 
related to accuracy. CTH is an artificial viscosity code, and numeri- 
cally widens shock discontinuities. Generally speaking, at very long 
ranges and times the numerical dissipation in the code will cause 
inaccuracies in the almost linear wave propagation that is expected 
to occur. In particular, the numerical dissipation will attenuate 
stress waves at great propagation distances. In addition, resolution 
of these low amplitude waves as they propagate through regions of 
large hydrostatic compression in the Earth’s interior must be main- 
tained. While we have still studied the overall response of the 
antipodal region to the KT impact, the resulting amplitudes and 
strain due to the long-term wave interactions can be in error. None- 
theless, this approach has also been pursued by Watts et al. [ 1 0], and 
has led to useful investigation of the response of the antipodal 
region. Our numerical resolution is significantly finer than that of 
Watts et al. [10]. Ultimately, we will improve the accuracy of our 
numerical modeling by linking earlier time CTH results to the more 
accurate seismological analyses described earlier. 
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Tsunamis are impulse-generated waves and can be generated by 
any event that displaces a significant amount of water in a short 
period of time, on the scale of seconds to minutes. The most typical 
source of tsunamis are subduction zone earthquakes, when deforma- 
tion associated with low-angle thrusting displaces water during the 
event. Tsunamis are also generated by volcanic events and subaque- 
ous landslides; bolide impacts in water will produce the largest 
tsunamis in Earth history. Tsunamis are wave trains, rather than 
single waves, and they may be locally damped, focused, amplified, 
and reflected. 

Our knowledge of typical (seismogenic) tsunamis and their 
deposits has increased significantly since Snowbird II. These tsuna- 
mis typically generate near field run-ups on the scale of meters, to 
tens of meters in rare cases. Recent tsunamis in Nicaragua (1992), 
Indonesia (1992), and Japan (1983, 1993) have allowed workers to 
generate detailed information about the nature of these tsunamis, 
their erosional an_d other destructive effects, and their deposits. 
Paleotsunami deposits from subduction-zone earthquakes have been 
documented from Chile (1960), Alaska (1964), the Pacific North- 
west U.S. (late 1600s), and elsewhere. I will present a brief sum- 
mary of the Nicaraguan, Chilean, and Pacific Northwest examples. 

Tsunamis generated by bolide impacts (“megatsunamis”) are at 
the other end of the tsunami scale, potentially generating tsunamis 
several orders of magnitude more energetic than seismogenic cases. 
Our understanding of these kinds of events is much more theoretical 
because there are no historic cases. We can expect that an impact in 
deep water could propagate a wave equal to the water depth, and 
thus megatsunamis may affect the floor of the ocean at depths much 
greater than seismogenic tsunamis, which only have an effect in the 
coastal zone. Certainly shelf environments (depths of -30,200 m) 



